Nonperturbative calculation of Born-Infeld effects on the Schrodinger spectrum of the 

hydrogen atom 
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We present the first nonperturbative calculations of the nonrelativistic hydrogen spectrum as pre- 
dicted by first-quantized nonlinear Maxwell-Born-Infeld electrodynamics with point charges. Judged 
against empirical data our results significantly restrict the range of viable values of the new electro- 
magnetic constant /3 introduced by Born. We assess Born's own proposal for the value of j3. 
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In the twenty years since its rediscovery in the 26- 
dimensional bosonic string theory study by Fradkin and 
Tseytlin 0, the nonlinear electromagnetic field theory 
proposed by Born and Infeld has been experiencing 
an astonishing renaissance. Recent surveys are 0, El- 
Most investigations since have been conducted from 
the perspective of the high energy community and involve 
higher-dimensional versions of the Born-Infeld theory (as 
in III LI ) and/or non-commutative analogs of it (as in jfj). 
Inevitably this has rekindled the interest in the original 
four-dimensional theory, the subject of this letter. 

We recall that Born's agenda was to rid (early) 
QED from its ultraviolet divergencies by quantizing self- 
regularizing nonlinear classical field equations. It was 
noted already in that the nonlinear Maxwell-Born- 
Infeld field equations In Born proposed a simpler 
nonlinearity than in ; yet in the electrostatic limit both 
field theories coincide, do not lead to the infinite self- 
energy problems of a point charge which occur with the 
linear Maxwell-Lorentz field equations, but the nonlin- 
earity made it difficult to proceed. With the spectac- 
ular quantitative successes of renormalized QED since 
the late 1940s, Born's original motivation became obso- 
lete; or so it would seem. However, as emphasized by 
Weinberg 0, more than half a century later standard 
QED is still in need of extrinsic mathematical regularizers 
to overcome the infinite self-energy problems of a point 
charge that have been inherited, in a sense, from the clas- 
sical Maxwell-Lorentz electrodynamics. In view of this, 
Born's suggestion ,7| to pursue some intrinsically self- 
regularizing nonlinear electromagnetic field theory reads 
as contemporary as it did in the 1930s; the rediscovery 
of Born-Infeld type Lagrangians in string theory, which 
could hardly have been foreseen by its founders, makes 
Born's suggestion all the more prophetic. 

The avoidance of infinite self-energies, as well as some 
other conceptual items are greatly to the theory's 
credit but surprisingly little is known about the empirical 
validity of the Born-Infeld theory. While the theory does 
not seem to have problems at the classical level [lj, [13 
it remains to be seen whether it will live up to its expec- 
tations at the quantum level. 

In this vein, a very natural question to ask is the follow- 
ing: What (detectable) effects does a hypothetical Born- 
Infeld nonlinearity of the electromagnetic fields have on 



the atomic spectra? This question should have been an- 
swered long ago. It was not, presumably because the 
nonlinearity of the field equations causes "difficulties [...] 
with the passage to the quantum theory, which appear to 
be insoluble with present methods of quantization" [l3j | 
(p. 32), and by 1969 "[t]he adaption [...] to the principles 
of quantum theory and the introduction of the spin ha[d] 
[...] met with no success" (p-375). As long as this 
situation prevails, one has to settle for quantum mechan- 
ical computations of spectral data in which Born-Infeld 
effects can be incorporated through the classical fields. 

Unfortunately, because the complicated nonlinearity 
of the field equations has stood in the way of finding 
relevant solutions with two or more point charges, all 
previous attempts to compute such quantum mechanical 
spectra 0, Jig, 0] have been foiled. In 0] the elec- 
tron is treated as a test particle in the known (see 0) 
Maxwell-Born-Infeld field of a point nucleus to compute 
hydrogen-like Schrodinger spectra to first order in pertur- 
bation theory; however, as we will see in this article, test 
particle theory is misleading for Born-Infeld equations. 
In 0] and 0], which have become standard references 
(see the introductions in 0] and (lflj). Dirac spectra are 
computed without recourse to test particle theory (al- 
beit with other approximations which are not of concern 
here) , defining the interaction energy as difference of the 
electrostatic field energy integrals for the bound versus 
the free configurations. However, the authors of and 
[l7j . who use Coulomb's solution Dq of the displace- 
ment field equation V • D — Anp with a charge density 
p comprising a single spectral electron and a spherically 
symmetric nucleus of charge z and a Thomas- Fermi cloud 
for the remaining z — 1 electrons, fail to realize that the 
nonlinear Born-Infeld law for the electromagnetic vac- 
uum maps this Coulomb field Dq into an electric field 
■Efgrs = Fbi{Dq) which is not identically curl- free [2oT |: 
more precisely, V x Tbi{Dc) 7^ almost everywhere, 
invalidating the spectral results of 0] and [l7j|. 

Recently, a consistent first quantization of the non- 
linear Maxwell-Born-Infeld field equations with point 
charges was achieved using the electromagnetic poten- 
tials [2l|| . Moreover, an explicit integral formula for 
the electron's electrostatic potential in certain proton- 
electron configurations (treated as point charges) was de- 
rived; this integral formula is readily extended to nuclear 
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charges z > 1 (see below). Thus the stage has been set 
for a systematic investigation of the simplest atomic and 
ionic spectra, the hydrogen- like ones. 

In order to keep technical matters as simple as possi- 
ble, here we only address the non-relativistic Schrodinger 
equation of a spinless electron bound to an infinitely mas- 
sive point nucleus. We plan to deal with the fine details 
contributed by relativity, spin, and the finite mass and 
size of the nucleus elsewhere. Furthermore, detailed eval- 
uations of the interactions and the eigenvalues are carried 
out only for the hydrogen atom (z = 1); the details of 
hydrogen-like interactions and ionic spectra for nuclear 
charges z > 1 are beyond the scope of this letter. 

In units of h for both action and magnitude of angular 
momentum, elementary charge e for charge, electron rest 
mass m e for mass, speed of light c for velocity, and Comp- 
ton wave length of the electron Ac = h/m e c for both 
length and time, a hydrogen-like spectrum is determined 
by the following dimensionless stationary Schrodinger 
equation on the electron's configuration space [2l|. 



(1) 



where s e is the electron's generic configuration space co- 
ordinate and the subscript e on Vg indicates differen- 
tiation with respect to s e . The fine structure constant 
a = e 2 /hc 1/137.036 is the dimensionless electro- 
magnetic coupling constant for the dimensionless total 
electrostatic potential (ftp defined below. The positive 
parameter (3 is Bom's electromagnetic vacuum constant 
("aether constant" for short) |22j, which enters through 
the Born-Infeld aether law, relating E (and H) with D 
(and B). Born argued that (3 — (3b with 



f3 B ~ 1.2361a. 



(2) 



Our spectral results allow us to assess the viability of (J2J) 
and Born's reasoning for it. 

The total electrostatic potential <j>p{s) at the actual 
space point s is determined by the electrostatic Maxwell- 
Born- Infeld equation V • 3~bi(— V<^g) = inp with p con- 
sisting of one positive and one negative point charge with 
values z and —1 at generic positions s„ and s e , respec- 
tively; explicitly, 



4tt (z5 Sn (s) - 5 Se (s)) , (3) 



with the asymptotic condition that 4>fi{ s ) ~ * for 
|s| — * oo. The solution of © depends on s as vari- 
able and on s n and s e as parameters; we sometimes em- 
phasize this by writing <f>i3(s\s n , s e ). While no explicit 
formula for <pp(s\s n , s e ) is known, Jpl reveals that we 
need to know only ^(s e |s„, s e ). Fortunately, although 
V x ^BiiDcis)) 7^ for almost every s in space, we 
do have V x !Fbi(Dc( s )) — for all s on the straight 
line through the point charges (this result generalizes to 



the vanishing of V x JF BI (Dc(s)) on the straight line 
through the respective centers of any two spherically 
symmetric charge distributions). Hence, an electrostatic 
potential function </>(s) solving for space points s 
on that line can be computed through the line integral 
(f)(s) = Tbi{Dc ( s ')) • ds' . Assuming D = Dq in 
leading order in (3, on this line we can approximately set 
(f>(s e ) = 4>j3(s e ). For z > 1 the integral is formidable, but 
when the nucleus is a proton (z — 1 and s„ = s p ), it can 
be recast into the more managable form [21| 



3 e \ a p , J 



f'(y) 



dx , (4) 



where r = \s p — s e |, xy = (3/r, and /' is the derivative of 



y 2 - y 



(5) 



For the remainder of this letter, z = 1. 

A look at the integral makes it plain that 
(3(t>p(s e \s pi s e ) depends on s p , s e and (3 only through the 
combination \s p — s e \/(3; hence, (3<j)p{s e \s p ,s e ) —: W(r/(3) 
is a function of r/(3. And while W does not seem to be 
expressible in terms of known functions, 10} lends itself 
readily to an analysis when the electron is far from, re- 
spectively near the proton. Note that "far" and "near" 
are relative to (3. 

If the electron is far from the nucleus, i.e., if r > 2v / 2/3, 
then W(r/P) can be expanded in an asymptotic series in 
powers of (3/r (asymptotically exact to four orders as 
r — > oo), thus 



W(r/f3)=Y. b x W/r) k + o(((3/r 



(6) 



fc=0 



with b = -\B (i, i), h = 1, b 2 = \B (|, |), and b 3 = 
2, where B( . , . ) is Euler's Beta function. Formula © 
reveals three important results. First, when r — > oo the 
electric potential at the location of the electron, 4>p(s e ), 
converges to the finite electron self-potential in Born- 
Infeld theory, defined by setting s = s e in Born's solution 



(s|s e 



d.z: 



A/0 



(7) 



for the electrostatic potential at s generated by a single 
(negative) unit point charge at s e 0, (NB: solves 
© when z = 0). We recall 0j: there is no short dis- 
tance Coulomb singularity of the single particle potential 
in Born-Infeld theory. Second, to leading order for large 
separation of electron and proton, the potential 4>p(s e ) 
varies with r reciprocally, i.e., we recover Coulomb's 
law for the pair potential from the nonlinear field equa- 
tion J3J. Third, there are higher order corrections to 
Coulomb's law. Indeed, when the electron is near the 
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nucleus, deviations from Coulomb's law become signifi- 
cant. More precisely, for r < 2^/2(3, the function W(r/(3) 
can be expanded into a Taylor series in powers of r/f3, 

oo 

W(r/P)=Y^*k{r/P) ik+1 (8) 

fc=0 

with explicitly computable expansion coefficients au ■ The 
first four of them read as follows: oq = —1/2, a\ — 
3/40 - 3tt/138, a 2 = -29/672 + 225tt/16384, and a 3 = 
1667/54912- 20265tt/2097152. Note in particular that 
W(0) = 0: there is no short distance Coulomb singularity 
of the pair potential in Born-Infeld theory [2lJ • 

Our discussion of (f>p supplies all the information we 
need to solve the Schrodinger equation (JIJ. To facili- 
tate the comparison with the familiar Schrodinger equa- 
tion for the Coulomb interaction, we write the eigenval- 
ues as E = % \B (j, jj + e and the total potential as 

Mse) = -$\B (|, \) + where Z(r/(3) is the ef- 

fective Coulomb charge of the proton "seen" from a dis- 
tance r. The self-potential terms on left and right hand 
side of JIJ then cancel out, leaving us to solve 

-IV^( Se )-aSrM^( Se ) = ^(s e ). (9) 

The Born-Infeld Schrodinger and the Coulomb 
Schrodinger potentials are compared in Fig. 1. 

Our first important spectral result states that the 
Coulomb limit (3 — > of eq. @ exists. In this case, 
Z(r/f3) — > 1 for all r > 0, which follows from 21]. 
Hence, in the limit (3 — > the spectrum of © converges 
to the familiar Rydberg law, i.e., 

*nl m = 71 = 1,2,..., (10) 

where n = 1,2, 3, ... and i = 0, 1, n — 1 and m = 
— £, ...,0, ...,£ are the usual main, secondary, and mag- 
netic quantum numbers. As is well known, t and m do 
not contribute to the energy eigenvalues m , so that 
n 2 of them coincide; we recall that this high degeneracy 
is due to the 0(4) invariance of © when /3 = 0. 

For all < (3 < oo, the 0(4) invariance is broken and 
the energy eigenvalues in general display only the 



FIG. 2: e /a 2 vs. a/3. 

2£+ 1-fold degeneracy corresponding to the manifest 0(3) 
invariance of JSJ; i.e., does not depend on the quan- 

tum number m. The 0(3) symmetry allows us to treat 
© by the usual separation of variables. Shifting the ori- 
gin of space to s p , the electron-proton distance r becomes 
the radial variable of standard spherical coordinates 
r, ip. In these coordinates the eigen-wavefunctions 
take the form V>S, m 0e) = R^Ji^Yf 1 ^,?), where the 

Y £ m (t?, </?) are spherical harmonics, and the R„l(r) satisfy 
the Sturm-Liouvillc problem 

(r 2 R r y-ie(£+l)-2arZ(r/(3)-2er 2 ]R = (11) 

for Jq 00 r 2 R 2 {r)dr < oo. We solved this radial problem 
by standard shooting technique, using MAPLE's Rungc- 
Kutta-Fehlberg45 method. 

It is instructive to discuss first the dependence of the 
ground state energy e ({3) = s[^ on {3. In Fig. 2 we 
display our numerically computed values of £${(3) for a 
selection of (3 values vs. (3, together with semi-explicit 
upper and lower bounds on £o (J3) , computed analytically 
except for numerical quadratures. 

There are several remarkable features visible in Fig. 2. 
First of all, there are two values of (3 at which the ground 
state energy Eq({3) coincides with the familiar Coulomb 
value —a 2 /2. Thus there are two regimes where the 
Born-Infeld theory yields a binding ( or ionization ) en- 
er 9U (~ £ o((3)) compatible with the empirical data: a per- 
turbative one near (3 = (the obvious one, as already 
noted) and a highly non-obvious and non-perturbative 
regime near (3 ~ 1.83297/a. Too far away from these 
two values, the binding energy would be either unrealis- 
tically large or small. Between these two (3 values, the 
binding energy is enhanced compared to a 2 /2, reaching 
a maximum at about j3 ~ 0.24774/a, while to the right 
of (3 » 1.83297/a, the binding energy is diminished, con- 
verging to zero as (3 — > oo. We remark that our semi- 
explicit upper and lower bounds allow us to rigorously 
prove this nonmonotonic behavior of the binding energy. 

We emphasize that the nonmonotonic behavior of 
£o(J3) is a quite nontrivial result; in particular, there is 
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no hint of it when the electron is treated (in "first ap- 
proximation" ) as a test particle which "feels" only the 
proton's electrostatic potential computed from the Born- 
Infeld equations with a single point source, neglecting the 
electron's own feedback [15|; i.e., (j>p{s e ) = 4>p(s e \s p , s e ) 
in © is replaced by Born's Q solution for a positive 
point charge, </4 +) (s e |s p ) = ~^ B \s p \s e ). For r large, 
</>^(s e |s p ) ~ 1/r, too. But since \J\ + x 4 > x 2 , it fol- 
lows from that ^g + '(s e |s p ) < 1/r, so test particle 
theory predicts a diminished binding energy for all [3. 

For a judicious selection of (3 values we computed sev- 
eral higher eigenvalues. Of particular interest are Born's 
value (J2J and the value (3 ~ 1.83297/a where the binding 
energy coincides with the Coulomb value a 2 /2 at (3 — 0. 
Listed in the table below are the energies of the ground, 
the first excited s, and the first p states for these 8 val- 
ues, as well as the corresponding empirical data |23j. We 
display a/3 rather than f3, and —e/a 2 rather than e; also, 
we suppress the magnetic quantum number m: 



a(3 
0.0000 
6.6 xlO" 5 
1.83297 

empirical 



sfl/a 2 
0.50000 
0.50016 
0.50000 
0.49973 



-08) 
-2,0 



/a 2 
0.12500 
0.12502 
0.19766 
0.12493 



-45/« 2 

0.12500 
0.19101 
0.36737 
0.12493 



From the table and Fig. 2 we are able to delineate 
the physically viable range of (3 values. By inspect- 
ing the excited energies, we can immediately rule out 
(3 as 1.83297/a, and by continuity also its neighborhood. 
That leaves only the perturbative regime of sufficiently 
small (3. But how small is "sufficiently small"? In par- 
ticular, is (3b "small enough"? 

The second row in the table lists spectral data for (3 = 
(3b- Note that — e (/?B)/a 2 deviates from — e (0)/a 2 = 
1/2 (first row) by 1.6 x 10~ 4 , and from the empirical 
data even by 4.3 x 10~ 4 . Of course, — £o(0)/a 2 differs 
itself from the empirical data by 2.7 x 10~ 4 , but as is well 
known, after correcting — £o(0)/a 2 for the finite mass of 
the proton, the difference to the empirical data reduces 
to only —3.14 x 10~ 6 , and the agreement improves even 
more with relativistic corrections. It is therefore to be 
expected that even after correcting — e (/3B)/a 2 for the 
finite proton mass and relativistic effects, the difference 
to empirical data will remain at about 10~ 4 . 

Even more dramatic is the splitting of 0.066a 2 between 
the 2s and 2p energies computed with (J2J , which is a fac- 
tor 10 4 bigger than the 2p 3 / 2 — 2pi/2 "fine structure" 
(which is not even visible at the level of precision in our 
table). Hence, even from the spectrum of the simplest 
atom we conclude that Born's value is not physically 
viable! 

Pending verification of our results through a more re- 
fined treatment, viable values of (3, as far as spectral re- 
sults go, must be much smaller than a. We plan to study 



just how small (3 must be using a relativistic theory with 
spin. But would the elimination of J5J not be a bearer 
of bad tidings for the Born-Infeld theory? Not yet! Born 
did not use detectable energy differences to compute lf2"| 
but equated the static field energy of a point charge to 
the electron's empirical rest energy, which yields 



^/(Vi-^Iv^bI 2 -l)d 3 s = l 



(12) 



in our units of m c 2 = 1. The integral equals 
giving (J2J. Unless i|12|) can be tied to a dynamical con- 
cept, such as scattering of an electron, the elimination of 
(0) by our spectral results is not bad news for the Born- 
Infeld theory. But this clearly calls for a deeper inquiry. 
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